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ABSTRACT 



We explore the minimum distance from a host star for an exoplanet to be 
potentially habitable, in order to maximize future chances of finding other hab- 
itable worlds. We find that the inner edge of the Habitable Zone (HZ) for hot 
desert worlds is at 0.5 AU around a solar-like star (well within the orbit of Venus). 
The relative humidity is the key controlling factor in determining the inner edge 
distance because water vapor has a strong impact on the greenhouse warming 
of the atmosphere, yet too little water vapor will deactivate precipitation and 
enable CO2 to accumulate. We estimate that a relative humidity as low as 1% 
can be sufficient to maintain a liquid water cycle and wash out CO2 from the 
atmosphere. If the surface pressure is too low (~0.1 bar), the water loss timescale 
of the planet is too short to support life. If the surface pressure is too high (~100 
bars), we show using atmospheric circulation arguments, that the day-night side 
temperature difference on slow rotators and tidally locked planets is too small 
to enable an active water cycle. In contrast, the temperature difference on fast 
rotators with high surface pressure can be large enough to produce rain. In- 
termediate surface pressures (~1-10 bars) can provide suitable conditions for a 
water cycle independent of the planetary rotation period. We additionally find 
that the water loss timescale is infiuenced by the atmospheric CO2 level, because 
it indirectly infiuences the stratospheric water mixing ratio. If the CO2 mix- 
ing ratio of dry planets at the inner edge is smaller than 10~^, the water loss 
timescale is ~1 billion years, which may be too short for complex life to evolve. 
Hot desert worlds might be identified using transmission spectroscopy and atmo- 
spheric retrieval methods. A hot desert world has small amounts of greenhouse 
gases (water and CO2) in the stratosphere, and the background gas must be opti- 
cally inactive in the mid-infrared. The new inner edge limit is significantly closer 
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to the host star than previous estimates. Therefore, the number of potentially 
habitable confirmed/validated planets increases by a factor of 2-3. 

Subject headings: ... 
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Introduction 



One of the main driver of exoplanet sciences is to find habitable or inhabited 
exoplanets. By definition an exoplanet is habitable if liquid water, a crucial ingredient 
for life, is present on its surface, and if th e surface te mperature and pressure are such 



that complex organic molecules are stable (JLang 



19861). Habitability is inferred remotely 
by studying the atmosphere of exoplanets and retr ieving atmospheric physi cal conditions 



such as the temperature profile and composition (IBenneke fc Seager 



speci f ically the presence 



2005 



Kaltenegger et al 



of biosignature gases (e . g . 



2010 : 



Sterzik et al 



Pes Marais et al. 



2002 



20121 ) , or more 



Segura et al. 



2OI2I ). Surface properties might be retrievable if 



the atmosphere is optically thin in so-called window regions of the spectrum. For example, 
if the surface pressure is small a nd the atmosphere is not fully cover e d by clouds, the surfac e 



prop erties are retrievable (e.g.. 



2OI3I). 



Seager et al. 



2005 



Fujii et al. 



2011 : 



Hegde fc Kaltenegger 



We summarize highlights of previous work on the inner edge of the Habitable zone 
going from the smallest to the largest inner edge distance. Typically an Earth-like planet 
with an Earth-like atmosphere is studied, by changing one attribute at a time. We describe 
this attribute and the level of model complexity. The inner edge around a solar-like star is 



located at 0.77 AU 



:or an Earth-like p 



circulation models (Abe et al 



2OIII). 



anet with a limited water reservoir based on 3D global 



Vladilo et al 



( 120131 ) showed using an energy balance 



model and varying the surface pressure between 0.3 and 3 bars that the inner edge is located 
at 0.85 AU within the parameter range where their model is calibrated. By maximizing 
the cooling effect of clouds on an Earth-like planet, the inner edge could be located as 



close as 0.87 AU around the Sun based on ID models with clouds fiKitzmann et al. 



2OI0I). 



ID cloud-free ca 



with 10 bar N2 (JKasting et al. 



culations show that the inner edge is at 0.88 AU for an Earth-like planet 



19931 ). These calculations have been redone with updated 
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opacity data (JKopparapu et al.ll2013[ ). The new results show that the inner edge is 0.92 AU 



assuming an Earth-hke atmosphere on a super-Earth planec. 

Here, we find that the inner edge is much closer to the host star, when non-Earth-like 
atmosphere conditions are considered. Super-Earth and Earth-sized planets are expected to 
be diverse and an extension of diversity to atmospheres must be considered in habitability 
studies, if we want to maximize our chances of finding a habitable world. We must first 
identify the key parameters influencing habitability and then establish the broadest (but 
plausible) range where the parameters could reside. In this work, the extreme value is 
adopted for some parameters, and the full range is considered for others. We use a variety of 
methods to estimate habitability: we calculate ID radiative-convective temperature-pressure 
and mixing ratio profiles; we study the effect of atmospheric circulation on the climate using 
order of magnitude estimates; and we also calculate the water loss timescale (the length 
of time over which the planet loses its water reservoir). All these methods are combined 
to identify the atmospheric scenarios that minimize the inner edge distance, and to reject 
scenarios where the planet is not hospitable for life. 

The inner edge of the habitable zone is controlled mainly by the properties of the 
host star, the planet atmosphere's relative humidity, CO2 mixing ratio, the most abundant 
gas, and the surface pressure. The first order effect of the host star is related to the 
inverse square law. The stellar luminosity increase s with stellar mass, therefore the 



I993I ). The atmospheric 



habitable zone is more distant around massive stars (iKasting et al. 
relative humidity is controlled by the water cycle. The CO2 mixing ratio is affected by 
the carbon-silicate cycle as well as the water cycle (both cycles are discussed in detail in 
the next paragraph). The most abundant gas of close-in habitable exoplanets is expected 
to be optically inactive, e.g., not a greenhouse gas, and not H2 because H2 is a potent 



^The calculation with the 10 bar N2 atmosphere was not performed in the new study. 
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greenhouse gas due to collision induced absorption (jPierrehumbert fc Gaidosll201ll ). The 
surface pressure is determined by various processes including (but not limited to) accretion 
from the protoplanetary disk, outgassing, at mospheric escape, accretion /erosion during 



impacts. Other aspect s such as e.g., clouds 



2011 



Zsom et al. 



2012 



low surface pressures ( 



, e.g.,), tidal heating 



Joshi fc Haberle 



Kitzmann et al. 



Barnes et al. 



2010; 



Goldblatt fc Zahnle 



20121 ). and the surface albedo at 



2OI2I ) could also impact the inner edge distance. 



The water cycle and the carbon-silicate cycle are the two most important processes 
that can alter the atmospheric composition and the surface pressure, and so we discuss it in 
more detail. The water cycle provides a positive climate feedback. As the surface becomes 
warmer, the evaporation rate of the surface waters increase and water vapor builds up in 
the atmosphere. As the optical depth of the atmosphere increases with the water mixing 
ratio, the greenhouse effect becomes stronger promoting further warming and evaporation 



[see e.g. 



Goldblatt fc Watson 



2012 



and references therein). The carb on-silicate eye 
climate stabilizer feedback process operating on a 10® years timescale (jWalker et al 



e IS a 



198l|). 



CO2 is removed 



of silicate rocks ( lAbbot et al. 



Tom the atmosphe re by liquid water precipitation and surface weathering 



20121 ). The weathered rocks are transported to oceans where 
abiotic or biotic processes produce CaCOa that sediments to the sea floor. If plate tectonics 
is present on the exoplanet, CaCOs is subducted to the mantle and removed from the 
surface regions. CO2 reforms by metamorphism in the mantle and it is outgassed by 
volcanism that closes the carbon-silicate cycle. If the climate cools, liquid precipitation (a 
sink for atmospheric CO2) is reduced, therefore CO2 builds up in the atmosphere and warms 
the surface to a new equilibrium temperature. The width of the Habitable Zone (HZ) is 



affect ed by these two cycles for an Earth-like planet (JKasting et al. 
2018k 



1993 



Kopparapu et al. 



We do not model the outer edge of the Habitable Zone, because previous work suggests 
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that even free-floating planets could be habitable (IStevensod Il999l ). In other words, 
the outer edge of the Habitable Zone is at infinity. Such planets must have a thick H2 
atmosphere to reduce the outgoing thermal emission and allow for slow cooling. Habitable 
conditions can also be maintained for an extended period of time, if heat is generated in the 
interior of the planet by e.g., radioactive decay (Rogers&Seager, in prep.). 

The paper's outline is the following. In Sect. [2l we first describe the initial atmosphere 
profiles and how the convective tropospheric temperature-pressure profile is calculated 
(Sect. 12. ip . Then we introduce our climate code specifically developed for this study (Sect. 
12. 2p . We describe the infiuence of atmospheric parameters on the inner edge distance 
in Sect. 12.31 We develop atmospheric dynamics tools to assess the plausibility of best 
case scenarios in Sect. \2A\ and estimate the loss timescale of surface waters in Sect. 12.51 
Our results are outlined in Sect. |3l where we present our main findings and supporting 
arguments. In Sect. HI we discuss the formation and the observable signatures of hot desert 
worlds, and the model uncertainties. Finally we summarize our results and draw conclusions 
in Sect. [5l 

2. Methods 

Our aim is to find atmospheric configurations that provide habitable surface conditions 
as close as possible to the host star. We do this in three steps: first, we study how various 
atmospheric and planetary properties infiuence the location of the inner edge by calculating 
radiative-convective equilibrium temperature-pressure profiles (Sections 12. ![ 12. 2[ and l2.3p : 
second, we estimate the requirements for a water cycle because liquid water precipitation, 
surface weathering and volcanic outgassing regulates the atmospheric CO2 level (Sect. 12. 4p : 
and third, we estimate the water loss timescale to determine whether the surface water 
reservoir is stable on a multi-billion year timescale (Sect. 12. 5p . 



2.1. Initial atmospheric profiles and convection 

The T-P profile and volume mixing ratios are initialized with a convective zone that 
follows the moist adiabat from the surface to the tropopause and is prolonged with an 
isothermal stratosphere. The atmosphere is described by the surface conditions (pressure, 
temperature, and albedo), relative humidity, CO2 mixing ratio, background gas (typically 
N2), and planetary surface gravity. Unless otherwise stated, we use 60 layers per pressure 
decade. The top-of-atmosphere pressure is fixed at 10~^ bar because the atmosphere is 
optically thin in all spectral bands above this level. The surface pressure is varied between 
0.1 and 100 bars. First, we establish the mixing ratios of CO2 and the background gases on 
a dry pressure grid (without water vapor). We assume that these gases are well mixed in 
the atmosphere and their mixing ratios are constant as a function of pressure. Next, the 
moist adiabat is built from the surface upwards given t he dry surface pressur e (-Psurf.dry), 



temperature (Tgurf), and relative humidity ($) following 



Pierrehumbert 



torn : 



1) 



TH20 



dlnPdry Cp^bg 1 _L \^PzJh£. _L f L A L 1 r 

^ L ^P.bs \Rh2oT ^J Cp,igT\ J' 

where Rbg and RH2O are the specific gas constant of the background gas and water vapor, 
respectively; Cp^hg and Cp^H20 are the specific heat capacity of the background gas and water 
vapor, respectively; L is the latent heat of water vapor, and fH20 is the water vapor mass 
mixing ratio: 

JH20 p , l^j 

^Tlfyg -'dry 

where ruhg and mH20 are the molecular masses of the background gas and water vapor, 
respectively, and Psa.t(T) is the saturation vapor pressure at temperature T calculated 
according to the Smithsonian Meteorological Tables. The moist adiabat is followed until the 
temperature drops below the prescribed initial stratosphere temperature. The stratospheric 
water vapor mixing ratio is constant and it is determined by the tropopause temperature 
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and relative humidity. Finally, the total atmospheric pressure is calculated as the sum of 
the water vapor partial pressure and the dry pressure. 

2.2. Radiative-convective equilibrium profiles 

We want to determine the equilibrium planet-star separation for given surface 
temperature and pressure. Therefore, we develop a specific numerical approach to calculate 
the T-P profiles: we fix the surface temperature and determine the semi-major axis where 
the atmosphere is in radiative equilibrium. Our approach solves the inverse climate problem. 



because typically the exoplanet's semi- major axis is fixe d and the corresponding sur: 



temperature is calculated (see e.g., the EXO-P model of 



Kaltenegger fc Sasselov 



2OIOI). 



ace 



The equilibrium T-P profile is determined by iterating the following steps: 

• We calculate the fluxes propagating through each atmospheric layer in several 
wavelength bands, we treat the thermal emission of the exoplanet/atmosphere and 
the stellar irradiation separately, and evaluate the stellar irradiation fluxes at an 
arbitrary distance. 

• We determine the semi-major axis where the top-of-atmosphere fluxes are in 
equilibrium using the Bond albedo of the exoplanet and the outgoing long wave 
radiation. 

• The stellar irradiation fluxes are rescaled such that the incoming stellar flux 
corresponds to the value received at the new semi-major axis. 

• The stratosphere temperature profile is updated based on the net fluxes propagating 
through the layers. 

• The tropopause height is updated if necessary. 
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The iteration stops, when the T-P profile reaches equihbrium. These steps are described in 
the following. 



2.2.1. Radiative transfer 

Fluxes are calculated in a limited number of wavelength bands for computational 
efficiency. The bands range between 0.2 to 200 microns, and they are linearly spaced in 
log(A). Such spacing ensures that Planck functions of various temperatures are resolved 
with the same power. We use 400 wavelength bands. The sufficient number of bands were 
determined by convergence tests where the outgoing IR flux of a test atmosphere was 
monitored as the number of bands is gradually increased. 

We calculate the optical depth of each layer in each wavelength band given the mixing 
ratio proflles. Various absorption and scattering sources could contribute to the optical 
depth of a layer. Th ese are the greenhouse gas absorption (discussed in the next paragraph). 



Rayleigh scatte ring flHarvev et al. 



of O9 and Qq (ISander et al. 



20021 ). N2 flBorvsow 



1998 



Weber 



2003 



Sneep fc Ubachs 



2OO5I ). U y absorpt ion 



20061 ). c ollision induced absorptio n (CIA) of H2 (JBorysow 



20021), and CO2 f lGruszka fc 



Borvsow 



1997 



Pierrehumbert 



, and continuum absorption 



2011 



, Sect. 4.4.8). It is 



of water and CO2 (following the description of 

also possible to place cloud particles in the atmosphere. In that case, the absorption and 

scattering coefficients, and the asym metry parameter of spheri cal droplets are calculated 



and tabulated using the Mie code of 



Bohren k HuffmanI (119831 ). 



We use the l ine-by-line HITRAN database (IRothman et al. 



sums formalism (IWiscombe &: Evans 



1977 : 



Pierrehumbert 



20091 ) and the exponential 



20 111 ) to determine the band- 



averaged transmission through an atmospheric layer due to greenhouse gas absorption. 
The high resolution line-by-line absorption coefficient is calculated at a pressure and 
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temperature grid covering typical valu es encountered i n our atmosphere profiles. We use 
the semi-analytical Voigt line profile of [ 



Liu et al. 



using the Doppler and Lorentz line widths (jOlivero 



(|20011). a nd determine the Voigt line width 



19771 ). The histogram of absorption 



coefficients are pre-calculated and tabulated, and these tables are read and interpolated at 
run-time to the pressure and temperature values of each layer. 

Once the optical depth of the atmosphere is known, fluxes propagating through the 
layer interfaces are calculated. We use the two-stream approximation and adopt the 
quadrature method for th e stellar illuminat ion, and the hemispheric mean approximation 



for the thermal radiation (JToon et al 



19891 ). The thermal flux of the planetary surface and 
atmosphere are treated separately from stellar irradiation. The distinction of the two fluxes 
is important to determine the equilibrium distance of hot inner edge planets with Ts^ri ~ 
400-500 K orbiting low temperature M stars {Tcs ~ 2000-3000 K). The outgoing thermal 
emission of the planet and atmosphere might overlap with the reflected stellar light in the 
near-IR. However, it is necessary to determine the relative contribution of these fluxes 
determine the equilibrium distance of the planet (see Eq. |3]). 



2.2.2. Climate model 

The goal of the climate code is to iterate the temperature and water vapor profiles 
to a radiative-convective equilibrium based on the fluxes propagating through the layer 
interfaces. This is achieved in three stages. First, we determine the planet-star separation 
where the top-of-atmosphere incoming and outgoing fluxes are in equilibrium, and rescale 
the stellar fluxes using the new distance. Then we update the stratospheric temperature 
profile, and finally the tropopause height is determined. All steps are outlined below. 
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The distance where the top-of-atmosphere fluxes are balanced is calculated as 

OTE = ^{l-a)FiM, (3) 

where OTE is the top-of-atmosphere outgoing thermal emission, a is the Bond albedo 
of the planet, and Fin is the top-of-atmosphere incoming stellar irradiation. As the T-P 
and mixing ratio profiles are assumed constant within one iteration, the outgoing thermal 
emission is independent of distance a. The albedo (the ratio of top-of-atmosphere outgoing 
and incoming stellar fluxes) is also distance- invariant, and it is evaluated at an arbitrary 
distance typically where the incoming stellar flux equals the Solar constant. Therefore, 
the unknown in Eq. [3] is Fj„ and the distance a where Fj„ is received at the top of the 
atmosphere. Once the equilibrium incoming stellar flux is determined, the downwelling and 
upwelling stellar fluxes are renormalized. 

We update the stratospheric temperature in the next step. We adopt a numerically 

efficient and simple approach to update the T-P profile: the temperature of a layer is 

increased by rfT, if the net flux propagating through the layer {dFi - the difference between 

incoming and outgoing fluxes) is positive, and the temperature decreases otherwise. The 

temperature step {dT) is 0.5 K in our simulations. Our method differs from traditional 

approaches that calculate the derivative dTi/dt (change of temperature over time) in each 

atmospheric layer i: 

dTj ^ _^dF\ 

where g is the gravitational constant, Cpi is the total heat capacity of the layer, dFi is the 
net flux propagating through the layer, dPi is the pressure difference at the top and bottom 
of the layer. If the net flux is negative, the layer cools. The time step [dt) is adaptively 
determined by setting an upper limit on dTi (dTmax)- We find that although the traditional 
approach properly captures how the temperature profile evolves in time, it converges slowly 
if the tropopause pressure level is orders of magnitude larger than the top-of-atmosphere 
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pressure level. Furthermore, it is unimportant how the temperature profile evolves in time 
from an arbitrary initial state, as the goal is to determine an equilibrium profile. 

Finally, the tropopause height and the stratospheric water mixing ratio is updated 
in the climate code. The lapse rate between the tropopause layer and the layer above is 
calculated. The tropopause height is increased by one layer, if the lapse rate is steeper than 
the moist adiabat. The tropopause height is decreased by one layer, if the lapse rate is 
larger than K/km (i.e., the layer above the tropopause is warmer). As the tropopause 
temperature changes, so does the tropopause and stratospheric water mixing ratios. 

The steps outlined in the previous two sections, namely radiative transfer, deter- 
mining the equilibrium distance, stratosphere temperature, tropopause, and water mixing 
ratio updates, are iteratively performed until the temperature profile relaxes to equilibrium. 
Our convergence criteria measures the temperature change of each layer for n iterations. 
Convergence is reached, if the temperature change is less than 2 dT in all of the layers 
during n successive iterations, where n is typically 10. Model validation and tests are 
described in App. |Al 

2.3. The ranges of model parameters, and their effect on the inner edge 

We provide a qualitative description on how the various atmospheric parameters 
influence the outgoing thermal flux, the albedo of the planet, and thus the inner edge 
distance in this section. In general, increasing the outgoing thermal emission and/or the 
albedo of the planet moves the inner edge closer to the star because more incoming stellar 
flux is necessary to achieve radiative-convective equilibrium. The most important effects of 
these parameters is described in detail in Sect. [3l while other parameters are fixed during 
the simulations. 
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Stellar type We perform simulations using various s tellar types ranging fron a M6V to 



FOV, or from 3000 K to 7200 K effective temperatures (ICastelli fc Kurucz 



2004h . The first 



order effect of stellar type is related to the luminosity of main sequence stars. An exoplanet 
has to be much closer to an M dwarfs than to an F star to receive the same total flux at the 
top of its atmosphere. However, the spectral energy distribution of stars have higher order 



effects on the inner edge distance. These effects become app arent, if the top-q 



atmosphere 



20131 ). The 



fluxes are compared at the inner edge around various stars (iKopparapu et al.l 
spectral energy distributions (SEDs) of low mass stars peak at longer wavelengths where 
Rayleigh scattering is less efficient, thus the planetary albedo is reduced. Furthermore, 
a larger fraction of the incoming stellar radiation is directly absorbed by the atmosphere 
around low mass stars because the greenhouse gases such as H2O and CO2 absorb in the 
near JR. To a smaller extent, the UV flux of the host star (especially the extreme UV) 
influences the composition of the stratosphere, which in turn can slightly influence the 



surface climate (see Fig. 4a of 



Rugheimer et al. 



2OI2I). 



Surface gravity The effects of surface gravity are two-fold. If the surface pressure is 
fixed, the cumulative optical depth of the atmosphere is reduced at all wavelengths for a 
larger surface gravity because the atmospheric scale height and thus the column density of 
greenhouse gases are reduced. On one hand, this increases the outgoing thermal emission 
of the exoplanet and shifts the inner edge closer to the star. On the other hand, more 
stellar flux reaches the surface because less stellar flux is absorbed by the atmosphere, 
and the column density of Rayleigh scattering background gas is reduced. If the surface 
albedo is low, the Bond albedo of the planet is reduced and the inner edge is pushed further 
away from the star. If the surface is reflective, the inner edge is pushed closer to the star. 
Therefore the albedo-effect of surface gravity depends on the relative contributions of two 
effects to the Bond albedo: Rayleigh scattering, and surface reflection. The mass and 
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radius of an exoplan et also influences how much atmosphere can the planet retain (see e.g., 



Lammer et al. 



20091 ) , and what is the volatile content of the mantle that can potentially be 



outgassed. 

We use values between g = 5 m/s^ for a planet somewhat larger than Mars, 10 m/s^ 
for an Earth-like planet, and 25 m/s^, which is a typical surface gravity of a 10 Mgarth 
super-Earth. 

Surface pressure The effect of surface pressure on the inner edge distance is two-fold. 
Larger surface pressure increases the albedo of the planet due to Rayleigh scattering. 
However, an increasing surface pressure also raises the column density of the greenhouse 
gases and thus the cumulative optical depth of the atmosphere for a fixed the mixing 
ratio. Therefore, the effect of surface pressure on the inner edge distance cannot be readily 
estimated, thus we treat it as a free parameter. We consider planets with 0.1, 1, 10, and 
100 bars of surface pressure. 



CO2 mixing ratio The lower the CO2 mixing ratio is, the closer the inner edge of the 
HZ is to the star. CO2 in combination with water vapor is a very effective greenhouse gas 
because these gases have non-overlapping absorption bands in the infrared. If CO2 was 
removed from the atmo sphere of Earth , the g lobal average temperature would be below the 



freezing point of water (jPierrehumbert 



2OIII ). Therefore it is crucial to calculate how the 



mixing ratio of CO2 influences the inner edge distance. 

The atmospheric CO2 level is determined by the balance between the sources and 
sinks of CO2. The abiotic source of CO2 is outgassing from the planetary interior. CO2 is 
removed from the a tmosphere by liquid water precipitation followed by silicate weathering 



( jWalker et al. 



19811 ). Atmospheric CO2 is converted into carbonates during the weathering 
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process. The carbonates can either dissolve in water as they sediment at deep ocean 
regions, or accumulate in solid form at shallow water regions. The role of plate tectonics in 
the carbon-silicate cycle is to recycle carbon. Plate tectonics removes the carbonates from 
the shallow water regions by subduction and CO2 reforms at the high temperatures of the 
mantle. Later CO2 is outgassed by volcanism that closes the carbon-silicate cycle. 

If an exoplanet has no plate-tectonics but volcanism and surface weathering, carbonates 
can accumulate in the shallow water regions until the mantle becomes devoid of CO2. 
As the source of CO2 at this point is depleted, liquid precipitation in combination with 
surface weathering can remove all CO2 from the atmosphere. The timesc ale of complete 



atmo spheric CO2 removal on Earth is estimated to be 400 million years (JKasting et al. 
I993I ). The timescale over which the mantle becomes carbon- free depends on the CO2 
outgassing rate and the carbon content of the in terior. By yarious estimates, the mantle 



of Earth holds 6 x 10 ' - 4 x 10 ° kg of carbon ( iJavoy et al 



Coltice et al 



1982 



Sleep fc Zahnle 



2004). The outgassing rate of CO2 is 3.3 x 10" kg/year (IWalker et al 



2001 



19811 ) 



Given these estimates, the carbon depletion timescale of Earth's mantle could be anywhere 
between 10^ and 10^ years, which is at least an order of magnitude shorter than the 
main-sequence life of the Sun. Therefore an exoplanet with an active hydrological cycle but 
no plate tectonics could be C02-free, unless CO2 sources other than volcanism exist on the 
planet. 

The presence of plate tectonics on exoplanets is difficult to asses both observationally 
or theoretically. Our current understanding is that it is less likely to have plate tectonics 
on super-Earth exoplanets. However, it cannot be ruled out due to large uncertainties in 
modeling and the lack of experimental d ata on core /mantle mate rials at pressures higher 
than experienced at the center of Earth f.Stamenkovic et al.ii201l[ l. 



As both the sinks and the sources of CO2 are highly uncertain on exoplanets, it is not 



-17- 

possible to self-consistently calculate the atmospheric CO2 level on an exoplanet different 
from Earth. Therefore, the CO2 mixing ratio is a free parameter in our study. The mixing 
ratio of CO2 is varied between 10^^ and 10^^. Large mixing ratio is expected, if plate 
tectonics is present on the exoplanet, because there is a stable source of CO2. However, 
the sink of CO2, liquid water precipitation is expected to be less efficient than on Earth 
because of the low relative humidities on hot desert worlds. The atmospheric CO2 level 
is expected to be small, if plate tectonics does not operate on the exoplanet, because the 
mantle becomes depleted in carbon on a billion year timescale, and liquid precipitation in 
combination with surface weathering removes CO2 from the atmosphere. 

Relative humidity The relative humidity is a key parameter in our study because it has 
a strong influence on the inner edge distance. Unfortunately it is exceedingly difficult to 
model the relative humidity profile of an exoplanet because it is determined by complex 3D 
processes such as atmospheric circulation, cloud formation, precipitation. The water cycle 
is also influenced by the spatial distribution of liquid water on the planetary surface. We 
do not attempt to model these effects. However, we highlight some crucial aspects of the 
water cycle, and provide qualitative arguments on the requirements for low tropospheric 
relative humidity. 

Strong temperature differences are necessary to maintain low relative humidity in the 
troposphere. The source of water vapor in the atmosphere is evaporation from a surface 
liquid reservoir, the sink of water vapor is condensation followed by precipitation. If 
precipitation does not reach the surface, the atmosphere eventually becomes fully saturated 
even if the fraction of surface liquids is small. Condensation and precipitation is initiated 
by cooling atmospheric parcels to saturation. Therefore, low relative humidity can be 
maintained, if strong temperature differences exist in the troposphere either vertically or 
horizontally. We argue that horizontal temperature differences are more important, because 
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if an air parcel is advected horizontally to cooler regions (night side or poles), rain drops 
surely reach the surface. However, if a rising air parcel e.g., in the Hadley circulation 
becomes saturated and rain drops form, it is not guaranteed that those drops reach the 
surface. The drops could evaporate as they fall through gradually warmer and drier regions 
of the troposphere (i.e., virga). 

We show in Sect. 12.41 that the lower limit of plausible relative humidity is $ ~ 1%. The 
lower limit is set by the required presence of liquid precipitation on the planet to wash out 
CO2 from the atmosphere. If the relative humidity is $ < 1%, the temperature difference 
necessary to cause precipitation is too large, and precipitation occurs in the form of snow. 
Solid precipitation above the freezing point of CO2 is not able to remove CO2 from the 
atmosphere, thus the CO2 level and the surface temperature would rise on such a planet. 
We find that the temperature difference is safely in the liquid water precipitation range for 
all atmosphere types considered in this work, if the relative humidity is $ = 1%. 

Surface temperature The surface temperature at the inner edge of the Habitable Zone 
is determined by the surface pressure and the relative humidity. Traditionally the HZ is 
defined by the presence of liquid water on the surface. Therefore the inner edge temperature 
is given by the boiling point of water: 

T3,rf = max(r : P,,rf + ^^sat (T) > Psat (T) ) , (5) 

where Psurf is the dry surface pressure, ^Psa.t(T) is the partial pressure of water on the 
surface. The inner edge surface temperature of Earth is 410 K for a dry surface pressure of 
1 bar, and the surface relative humidity is 70%. If Earth had 1% relative humidity and all 
other properties kept constant, the inner edge surface temperature would be lower, 373 K 
because the partial pressure of water is reduced. 

The maximum surface temperature is limited by the chemical stability of complex 
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organic molecules. If only Eq. [5] is used to calculate the inner edge surface temperature, its 
allowed range is between 273 K and 647 K given by the triple and critical points of water, 
respectively. However, can life exist at temperatures as high as 647 K? On Ear th, organisms 



grow at temperatures up to 395 K (JKashefi fc Lovley 



20031 : 



Takai et al. 



20081). However, 



the DNA an d amino acids become chemically unstable only at temperatures above 500 K 



dLang 



19861 ). Therefore we limit the surface temperature to a maximum of 500 K. We note 



that this is a globally averaged surface temperature because we use a ID climate model. 



As discussed in the previous paragraph and in Sect. 12. 4[ large horizontal temperature 
differences are a requirement to initiate precipitation. Therefore the surface temperature at 
a substantial part of the planet must remain well below 500 K (see Sect. 12. 4p . 



Surface albedo We assume that our exoplanets are hot and dry, thus most of the surface 
at low zenith angles is possibly cover ed by deserts or barr en rocky surfaces. The surface 



albedo of Earth is 0.15 (Chapter 5 of 



Pierrehumbert 



2011 



based on ERBE measurements). 



and the main contri butors are wat er, land, vegetation, and ice /s now surfaces. Th e albedo 

19831 ). the albedo of water is 0.05 (Clark et al. 



of the Sahara is 0.4 (iTetzlafi 



20071 ) We do 



not expect large water-covered areas on hot desert worlds, therefore we expect surface 
albedos between 0.05 and 0.4. Our nominal value is 0.2, but we varied the surface albedo 
between and 0.6. We find that although higher surface albedos of low surface pressure 
atmospheres have an impact on the inner edge distance, the surface reflection becomes less 
important at high surface pressures because most of the incoming stellar flux is scattered 
back to space by Rayleigh scattering, or absorbed directly by the atmosphere. At 100 bars, 
typically less than 1% of the stellar flux reaches the surface. If the surface pressure is low, 
the surface albedo influences the tropopause height because the convective zone is driven 
by the absorption of stellar light at the surface. 
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Nominal parameter values The parameter values of our nominal simulation is shown 
in Tabled! The most important effect of these parameters is shown in Sect. [31 We typically 
change one parameter at a time, while the other parameters are held fixed at the value 
given in Tabled! 

Table 1: Overview of parameter values in the nominal simulation. 



Parameter Value 

surface gravity gsurf = 25 m/s^ 

surface pressure -Psurf = 1 bar 

CO2 mixing ratio Xco2 — 10""^ 
relative humidity <1> = 1% 

surface temperature Tsurf = 370 K ^ 
surface albedo Osurf = 0.2 

background gas N2 



2.4. The interplay of atmospheric circulation, and the hydrological cycle 

As discussed previously, liquid water precipitation is crucial for two reasons: to keep 
the relative humidity of the atmosphere low, and to reduce the atmospheric CO2 levels. An 
air parcel has to be cooled down to saturation to initiate condensation and precipitation. 
Therefore, the goal of this section is to study atmospheric circulation processes that can 
cause strong temperature differences. We distinguish two limiting cases: slowly rotating 
and/or tidally locked planets, and fast rotators like Earth. The distinction is necessary 
because different circulation regimes dominate the heat transfer depending on the rotation 



^The surface temperature is a function of relative humidity and surface pressure. 
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period. 

The dew point temperature (T^) necessary to initiate condensation and precipitation 
is given by the globally averaged temperature of the air above the surface and the relative 
humidity. As a point of reference, if the globally averaged surface temperature is 370 K, 
and the relative humidity is 1%, the parcel has to cool down to T^ = 280 K to reach 
saturation. The dew point temperature for the same relative humidity value is 337 K, if 
the surface temperature is 500 K. Dew point temperatures above 273 K ensure that liquid 
water precipitation is present on the planet to act as a sink for atmospheric CO2. 

Cloud formation by vertical mixing is not considered because precipitation from 
cumulus clouds might not reach the surface (i.e., virga). If condensation occurred e.g., on 
the upwelling part of the Hadley circulation around the equator, it is uncertain whether rain 
drops reach the surface. It is likely that the drops evaporate as they fall through gradually 
warmer and drier regions of the troposphere. 

2.4-1- Slow rotators and tidally locked planets 

The equator-to-pole temperature difference is typically small on planets with long 
periods (several tens to hundreds of days) because the Hadley cell extends all the way 

leat very efficient l y, the equator-to-pole 



1980 



Showman et al 



2OI1I) 



to the poles. As Hadley circulation transports 

temperature difference is expected to be small (iHeld fc Ho 

However, the day-night temperature difference can be significant, if the atmosphere 

efficiently cools on the night side. 

The speed of the zonal surface winds and the radiative cooling timescale determine 
the day-night temperature difference. Using the dew point temperature (T^), and the 
radiative cooling timescale of the atmosphere, the zonal advection timescale is estimated 
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( IShowman &: Guillot 



2OO2I Eq. 11): 



''"zonal 



Trarl l0£ 



T. 



surf 



T, 



T. 



surf 



where Trad is the radiative timescale of the atmosphere: 

AP c„ 

'rad 



^surf4aT3^rf' 



(6) 



(7) 



where AP is the pressure difference at the top and bottom of the layer above the surface, 
Cp is the heat capacity of the layer, and a is the Boltzmann constant. The zonal wind speed 
necessary to reach T^ is 

^zonal -'•'p/ ''"zonal; \^) 

where Rp is the planet radius. As the required temperature difference is large (~100 K), 
it is necessary to have a small wind speed to preserve the large temperature contrast - 
although temperature differences drive winds. Therefore, we rejects exoplanet scenarios for 
which the estimated zonal wind necessary to preserve the required temperature difference 
is too small and use t he mean surface wind speed of Earth as a reference value (7 m/s. 



Capps fc Zender 



20081). 



2.4-2. Fast rotators 



The day-night temperature difference on fast rotators is small, but the equator-to-pole 
temperature difference can be signific ant. The Hadley c ell cannot extend all the way to the 



poles on fast rotators, such as Earth ( JHeld fc Hou 



1980). Instead, heat is tran sported from 



the maximum latitude of the Hadley cell to the poles by baroclinic instability ( IStone 



19781). 



In energy balanc e models, baroclinic instability is approximated as a large scale diffusion 



process (see e.g. 



Held 



1999 



Spiegel et al. 



20081 ). We also take advantage of this formalism 



and calculate the diffusion coefficient necessary to create the equator-to-pole temperature 
difference of Tsurf — T^. 
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The diffusion equation of the energy balance model is (jShowman et al. 

ar(0) 



20111 ): 



dt 



V{cpDAT) + S{a, 0) - OTE, 



(9) 



where T(0) is the temperature as a function of latitude, D is the diffusion coefficient, and 
S{a,(j)) is the absorbed stellar flux as a function of latitude. An analytical solution exists 
for Eq. [9l if the albedo, D, Cp are constants, the OTE is fitted as A + BT, and if the stellar 
flux is para meterized a s a constant plus a term proportional to the Legendre polynomial 



P2(cos(0)) fiHeldlll999l ). The diffusion coefficient expressed from the analytical solution is: 



D 



J-suTl — J-d J P 

6c„ 



(10) 



The parameter B is determined by calculating radiative-convective equilibrium profiles 
with surface temperatures 1 K above and below the nominal Tgurf, and fitting the OTE. 
Our reference value for D is the diffusion coefficient that reproduces the zona lly averaged 



Suarez &: Held 



19791). If an 



equator-to-pole temperature profile of Earth (10^ ni^/s, 

atmospheric scenario requires D < 10^ m^/s to preserve the necessary temperature 

difference, we reject the scenario. 



2.5. The loss timescale of surface waters 

The water reservoir of the planet is gradually depleted because the stellar UV radiation 
dissociates the stratospheric water vapor, and hydrogen subsequently escapes to space. 
The goal of this section is to estimate the loss timescale of the exoplanet's liquid reservoir 
assuming diffusion-limited escape. If the water loss time is too short, the planet could loose 
its water reservoir before life develops. On the other hand, if the timescale is on the order of 
several billion years, it is safe to assume that the surface waters are stable and long-lasting. 

The most important reason to calculate radiative-convective T-P profiles is to reliably 
estimate the stratospheric water mixing ratio and through that the diffusion-limited escape 
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of H2O. The tropopause pressure, temperature and relative humidity set the stratospheric 
water mixing ratio. It is expected that both the tropopause pressure and temperature 
varies with atmospheric and stellar parameters. Often i t is assumed that the s tratosphere 



has a constant 200 K temperature (JKasting et al. 



1993; 



Kopparapu et al. 



20131). However, 



the water loss timescale might be inaccurate as the tropopause properties are not 
self-consistently calculated. We note that a constant 200 K stratosphere temperature is 
small enough not to significantly affect the OTE and thus the radiative equilibrium distance 
of the planet. 

We assume that our hot desert worlds have a 100 times smaller liquid water reservoir 
than Earth, thus the total mass of liquid water reservoir is 1.4 x 10^^ g, which amounts to 
A^HaO = 4.6 X 10"^*^ H2O molecules. The water loss timescale is calculated as 



(11) 



"" A top Tr ' 

^surf'^H20'^diff 

where v4surf is the surface area of the exoplanet, '".^^'o i^ ^^^ number density of H2O molecules 
at the to p of th e atmosphere, and Vdig is the diffusion velocity calculated according to 



Hu et al. 



(120121 ). If the top-of-atmosphere temperature is 200 K, and the water vapor 



mixing ratio is 3 x 10 ^ on a planet with 10 m/s^ surface gravity and one Earth radius, 
the water loss times c ale is 5.2 x 10^ years, in close agreement with the calculations of 



Kasting et al. 



(119931 ): 



Kopparapu et al 



mm . 



3. Results 

Our most significant finding is that exoplanets can be habitable much closer to the host 
star than previously estimated. For example, an exoplanet can be habitable as close as 0.5 
AU from a solar-like star given favorable atmospheric properties. Hot and dry desert worlds 
with 1% relative humidity, a broad range of CO2 mixing ratios and surface pressures, have 
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sufficiently long water loss timescales, such that surface hquid waters would continuously be 
present. The atmospheric circulation estimates also confirm this finding: the liquid water 
cycle can plausibly operate under a variety of atmosphere scenarios, thus CO2 is washed 
out of the atmosphere, and the relative humidity is kept low. In other words, a variety of 
atmospheres could keep a planet habitable at the inner edge. 

3.1. Relative humidity as the main controlling factor 

The inner edge of the Habitable Zone can be as close as 0.5 AU around a solar-like 
star. The relative humidity of the atmosphere is the most important factor controlling 
how close a planet can be to a star and still maintain surface liquid water. The relative 
humidity can be as low as 1%, a plausible minimum for the liquid water cycle to operate 
(as shown in Sec. 12. 4p . If the relative humidity is lower than 1%, the necessary temperature 
difference (from day to night side or from equator to pole) to initiate condensation would 
be too high and precipitation would occur in the form of snow instead of rain on the night 
side. As snow is not able to wash out CO2 from the atmosphere, CO2 would accumulate 
and a planet with volcanism would develop a strong greenhouse effect like the Venusian 
atmosphere. The habitable zone for planets of 1% and higher humidity levels is shown 
in Fig. [1] In the following we discuss the top of atmosphere radiative properties of the 
atmosphere: the outgoing thermal emission (OTE) and the albedo of the planet. 

The inner edge distance is close to the host star mainly because the outgoing thermal 
emission is high if the atmosphere is dry. In other words, the emitted thermal radiation 
escapes easily, if the greenhouse gas concentrations are low. The OTE does not depend on 
the stellar type (see Fig. |2^); it is mainly influenced by the pressure-temperature profile and 
the greenhouse gas concentrations (shown in Fig. |3]for the nominal atmosphere of Tabled]). 
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Fig. 1. — The inner edge of the Habitable Zone as a function of relative humidity. Other 
atmospheric parameters are fixed at the values given in Table [H The x axis represents the 
planet-star separation, the y axis shows the mass of the host s tar. The inner edge dist ance 



is small for low relative humidity. The Habitable Zone limits of 
indicated with black solid lines. 
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Fig. 2. — Left: The outgoing thermal flux as a function of relative humidity. The x axis 
shows the mass of the host star. Right: The Bond albedo of exoplanets as a function of 
relative humidity. The blue solid line indicates the surface albedo of 0.2. 
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Fig. 3. — The temperature and water mixing ratio profiles of the nominal atmosphere with 
parameters given in Table [TJ The color coding represents the T-P and water mixing ratio 
profiles of exoplanets orbiting various host stars. The effective temperatures of four selected 
host stars are indicated in the legend on the left side. The legend on the right side shows the 
equilibrium distance of the exoplanet around the host stars. The color coding is identical on 
all figures showing T-P and water mixing ratio profiles. 

The second factor influencing the inner edge distance is the albedo of the planet. The 
planetary albedo is str ongly influenced by the host star's spectral energy distribution, and 



the relative humidity (JMarley et al. 



19991 ). Low mass stars radiate at long wavelengths 
where absorption by greenhouse gases such as water and CO2 is present. Therefore, the 
incoming stellar light is taken up by the atmosphere before it reaches the surface. The larger 
the relative humidity, the more stellar light is absorbed (see Fig. 12b). On the other hand, 
massive stars radiate at short wavelengths where Rayleigh scattering dominates. Therefore 



the planetary albedo around massive stars is typically larger than the surface albedo. 
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3.2. CO2 mixing ratio influencing the water loss timescale 



The atmospheric CO2 mixing ratio influences the inner edge distance, which itself 
indirectly affects the water loss timescale. We discuss the effects of CO2 on habitability 
from two perspectives: in this work, changes in the CO2 level modifles the planet's radiative 
equilibrium distance (aka. the inner edge distance) because the surface temperature is flxed; 
the other perspective is to flx the planet's semi-major axis and investigate the influence of 
CO2 on the climate. We shortly discuss the second perspective at the end of this section. 
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Fig. 4. — Left: The temperature and water mixing ratio proflles for a CO2 mixing ratio of 
10~^ with 1 bar surface pressure, 1% relative humidity. Right: The temperature-pressure 
and water mixing ratio proflles for a CO2 mixing ratio of 10~^, other atmospheric parameters 
are unchanged. Higher levels of CO2 results in colder tropopause temperatures, and dryer 
stratospheres at the inner edge. The color coding is identical as in Fig. |3l 



The inner edge distance of planets with C02-rich atmospheres is large. The already 
knownreason is that CO2 is a potent greenhouse gas, thus the outgoing thermal emission 
decreases with higher CO2 levels. For example, the OTE is 800 W/m^ and 700 W/m^ for 
Xco2 = 10^^ and 10"^, respectively. However, the albedo remains largely unaffected. As a 
result, the inner edge distance is 0.48 AU, and 0.52 AU for CO2 mixing ratios of 10~^ and 
10~^ around a solar-like star. 
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Fig. 5. — The water loss timescale at the inner edge as a function of CO2 mixing ratio for 
atmospheres of 1% relative humidity and 1 bar surface pressure (for other parameters, see 
Table [1]). As the stratospheric water levels are higher for low CO2 mixing ratios (see Fig. 
H]), the water loss timescale becomes shorter. If the CO2 mixing ratio is less than 10~^, the 
water loss timescale becomes shorter than 10 billion years. 

The tropopause and stratosphere of exoplanets at the inner edge with low levels of 
atmospheric CO2 are warm (see Fig. |3]). This effect is explained by the variation of the 
inner edge distance with the CO2 mixing ratio. Close-in planets are strongly illuminated, 
therefore the tropopause and the stratosphere are warm. The tropopause temperature is 
240 K and 190 K for CO2 levels of 10"^ and 10"^ 

Close-in planets with low levels of CO2 might quickly lose their surface waters and 
become inhabitable. The water loss timescale as a function of CO2 mixing ratio is shown in 
Fig. |5l For a fixed relative humidity, the tropopause temperature determines the water-loss 
timescale (see Sect. 12. ip . As discussed in the previous paragraph, the troposphere is cold, 
if the atmosphere is C02-rich at the inner edge. Therefore, the atmospheric CO2 level 
has an indirect influence on the water loss timescale, because it regulates the tropopause 
temperature. 
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We note that for a fixed semi-major axis, the CO2 mixing ratio has an opposite 
(but small) effect on the water loss timescale. Radiative-convective models of Earth's 
atmosphere sho w that doubling the atmosph eric CO2 warms the surface but slightly cools 



the tropopause (ISchlesinger fc Mitchell 



19871 ). The result is that the stratospheric water 



mixing ratio reduces, and the water loss timescale is prolonged. 

3.3. Intermediate surface pressure necessary for habitability on dry planets 

The impact of surface pressure on habitability and the inner edge distance is complex. 

Low surface pressures (~0.1 bar) are a problem because the water loss timescale is 
short and complex life might not have enough time to evolve (see Fig. E]). The water loss 
timescale of a hot desert world is only 0.1 billion years if the surface pressure is 0.1 bar, this 
is because the troposphere is confined to low altitudes, and the tropopause temperature only 
slightly differs from the surface temperature (see Fig. [7^). Therefore the stratosphere is 
water rich and the water loss timescale is short. At large surface pressures, the troposphere 
is extended, the tropopause is cold (see Fig. Uh), and the water loss timescale is always in 
excess of 10 billion years (see Fig. |6]). 

High surface pressures (100 bars) on tidally locked planets or slow rotators are 
problematic for habitability. The liquid water cycle cannot plausibly operate because the 
day-night side temperature difference is too low to initiate condensation. The rotation 
period of the exoplanet is an important factor influencing the water cycle because it 
prescribes the dominant circulation mode that redistributes heat. On tidally locked planets 
or slow rotators, heat from the day side is transported by zonal flows to the night side 
(see Sect. I2.4.ip . If the surface pressure is 100 bars, the radiative cooling timescale of the 
atmosphere is long compared to the advection timescale. That is to say, an atmospheric 
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Fig. 7. — Left: The temperature and water mixing ratio profiles for an atmosphere with 0.1 
bar surface pressure. Other parameters are given in Table [1] Right: Atmospheric profiles 
for a hot desert world with 10 bar surface pressure. The tropopause is much colder in this 
case. The color coding is identical as in Fig. |3l 



parcel travels across the night side of the planet before it can cool down to the dew point 
temperature. To give a more qualitative argument, we calculate the surface zonal wind 
speed necessary to initiate condensation on the night side (Eq. |8]). The zonal wind speed is 
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on the order of 1000 m/s, 100 m/s, 10 m/s, and 1 m/s for Psurf = 0.1, 1, 10, and 100 bars, 
respectively. The inverse correlation between surface pressure and zonal wind speed is due 
to the long radiative cooling timescale at large pressures. The necessary zonal wind speed 
on exoplanets with 100 bars of surface pressure is too low, therefore, such atmospheres 
cannot plausibly cool down to the dew point temperature on the night side, and water 
condensation does not occur. On fast rotators like Earth, baroclinic instability transports 
heat from the equatorial regions to the poles (see Sect. 12.4.21) . We find that the baroclinic 
diffusion coefficient required to initiate condensation around the poles (Eq. [TOl) is on the 
order of 10^*^ - 10^^ m^/s for all surface pressures considered. These values are 4-5 orders of 
magnitude larger than the baroclinic diffusion coefficient on Earth. Therefore large equator 
to pole temperature differences and hquid precipitation are plausible on fast rotators at all 
surface pressures. The baroclinic diffusion coefficient is large because the OTE changes 
rapidly as a function of surface temperature in dry and hot atmospheres (parameter B in 
Eq. [10]). 
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Fig. 8. — The inner edge of the Habitable Zone as a function of surface pressure. The inner 
edge moves fu rther away from the host star with increasing surface pressures. The Habitable 



Zone limits of 



Kopparapu et al. 
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A habitable planet with the smallest semi-major axis should have a surface pressure 
around 1 bar for stars less massive than 1.2 M©, and a surface pressure of 10 bars around 
stars more massive than 1.2 Mq (Fig. [H]). The inner edge distance is a non-linear function 
of the surface pressure and the host star properties, because both the OTE and the albedo 
show complex features (Fig. [9]). The OTE is maximized if the surface pressure is 1 bar. 
However, water boils at low temperatures, if the pressure is reduced. Therefore the surface 
temperature and thus the OTE are small at low pressures. Although the boiling point and 
the surface temperature are high at large pressures, the optically thick atmosphere reduces 
the OTE at large pressures (see Fig. ^j^). The albedo of low pressure exoplanets is very 
close to the surface albedo. On the other hand, exoplanets with large surface pressures (100 
bars) simultaneously have the lowest albedo around low mass stars due to greenhouse gas 
absorption; and they also have the highest albedo around massive stars due to Rayleigh 
scattering (Fig. [Qjo). 
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Fig. 9. — Left: The outgoing thermal flux as a function of surface pressure. Right: The 
Bond albedo of exoplanets as a function of surface pressure. The x axis shows the mass of 
the host star. 



An exoplanet's surface pressure is actually difficult to ascertain whether it be a priori 
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from theory or from observations. The surface pressure of exoplanets is influenced by 
processes such as accretion from the protoplanetary disk, atmospheric escape, atmosphere 
accretion and/or erosion during impacts, and out gassing. In fact, the most challenging 
problem in proving that an exoplanet is habitable will be to observationally constrain its 
surface pressure. Even if the atmosphere is cloud-free, the atmosphere could be so optically 
thick that the surface remains hidden and unconstrained. 

3.4. Surface gravity weakly influences the inner edge 

The effect of the planet mass and radius on the inner edge distance is small compared 
to other factors discussed previously (see Fig. [10]). Surface gravity influences the pressure 
scale height of the atmosphere and thus the column density of gases for a fixed surface 
pressure. The atmosphere of a Mars-sized small planet is expanded and thus the column 
density of all gases is raised. This has implications for the OTE as well as the planetary 
albedo. Increasing the greenhouse gas column density reduces the OTE from 800 W/m^ 
to 700 W/m^ for planets with surface gravities of 25 m/s^ and 5 m/s^, respectively. The 
albedo of a Mars-sized small planet shows the largest deviation from the surface albedo 
(see Fig. [TT1) . The albedo around low mass stars is the lowest for small planets, because 
the greenhouse gas column density is large and the incoming stellar light is absorbed. At 
the same time, the albedo around massive stars is also the largest for small planets because 
Rayleigh scattering is more efficient due to the large column density of N2. 

Although the differences are small, the inner edge is slightly closer to the host star for 
super-Earth planets (see Fig. 8b). The distance around a solar-like star is 0.49 AU, 0.5 
AU, and 0.51 AU for surface gravities of 5, 10, and 25 m/s^. As the inner edge distance 
changes rather moderately with planet mass and size, it is a good approximation to use one 
common inner edge limit for all super-Earth and Earth-like planets. In the following we 
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Fig. 10. — The inner edge of the Habitable Zone as a function of surface gravity. The x axis 
is the planet-star separation, the y axis shows the mass of the host star. Surface gravity has 
a small effect on the inner edge distance for the atmospheric configuration adapted. The 
inner edge is s omewhat closer to the h ost star, if the surface gravity is small. The Habitable 



Zone limits of iKopparapu et al. 



( I2OI3I ) are shown with black solid lines. 



adopt the inner edge limit of the nominal atmosphere with properties outlined in Tabled]). 



3.5. Application to confirmed super-Earth exoplanets and candidates 

We show semi-major axes and stellar masses of the currently confirmed and candidate 
super-Earth exoplanets in rela t ion to our fundamental inner edge boundary and to the 



3Z limits of 



Kopparapu et al 



(|2013[ ). We use the Open Exoplanet Catalogue of 



Rein 



( 20121 ) as the source of stellar and exoplanetary data for confirmed exoplanets. The NASA 



Exoplanet Archive is used to access the Kepler exoplanet candidates of the Q1-Q8 data sets 



Burke et al. 



Torres et al. 



2013) 



. If missing, the stellar mass is calculated using the empirical relations of 



( I2OIOI ). We consider exoplanet (candidates) that have a mass (or msini) less 



than 10 Msarth, and/or a radius less than 2.5 -REarth- 
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Fig. 11. — The Bond albedo of exoplanets as a function of surface gravity (other parameters 
given in Tabled]). The outgoing thermal emission (not shown) is 700, 745, and 800 W/m^ 
for surface gravities of 5, 10, and 25 m/s^, respectively. As shown in Fig. [TOl surface gravity 
has a modest influence on the inner edge distance. 

There are nine radial velocity super-Earths, one confirmed transiting exoplanet, and 
one exoplanet discovered by microlensing that could potentially be habitable given suitable 
atmospheric conditions. The confirmed exoplanets are shown in Fig. [T2k . The potentially 
habitable radial v elocity planets with increasing host star mass are Gliese 581c, d, f, g 



( Udrv et al. 



20071. although the existence of Gliese 58 



( JAnglada-Escude et al 



HD 40307g fTuomi et al 



2012h. HD 85512b f lPepe et al. 



f an d g are debat ed), Gliese 667Cc 



201 ih. HD 20794d f|PeDe et al. 



20131 ). and tau Ceti d, e, and f (jTuomi et al 



potentially habitable exopla net is Kepler 22b ( 



exoplanet OGLE -05-390L b ( 



Beaulieu et a 



is hydrogen-rich (jPierrehumbert fc Gaidos 



Borucki et al 



201ll ) 



2OI2I ). The transiting 



2OI2I ). The microlensing 



2OO6I) could also be habitable, if its atmosphere 



20111 ) or if it has a sufficient internal heat source 



(Rogers et al, in prep.). 

We identify twenty potentially habitable Kepler candidates that have a radius below 
2.5 -REarth- Their KOI IDs in order of increasing stellar mass are K01686.01, K00854.01, 
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Fig. 12. — Left: The semi-major axes and stellar masses of confirmed super-Earth planets, 
and the fundamental inner edge of the Habitable Zone. Exoplanets discovered by vari- 
ous detection me t hods are included in the figure, as well as the Habitable Zone limits of 



Kopparapu et al. 



(120131 ) ■ and the positions of the four rocky Solar System planets. Right: 
The Kepler confirmed and candidate super-Earth planets. The number of potentially hab- 
itable exoplanets increases by a factor of 2-3 compared to previous estimates. However, 
follow-up observations aimed to characterize super-Earth atmospheres are necessary to de- 
termine whether these exoplanets are really habitable. 

K02418.01, K03010.01, K01361.01, K01298.02, K01174.01, K02882.01, K02770.01, 
K02762.01, K00518.03, K01871.01, K02834.01, K00701.03, K02469.01, K03036.01, 
K02474.01, K02311.01, K00172.02, and K01739.01. Although our minimum stellar mass is 
0.33 M0, Fig. [T2b indicates that potentially habitable exoplanets are quite possibly teeming 
around low mass stars. 
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4. Discussion 



4.1. The formation of desert worlds 



Where should we look for hot desert worlds? Is it possible to estimate their occurrence 
rate? What factors play a significant role in forming desert worlds? To answer these 



quest ions, we rely on N-body simulations of interacting planetary ernbryos and p 



[e.g. 



Morbidelli et al. 



2000 : 



Chambers 



2001 



Raymond et al. 



2004 . 



2006 



anetesimals 



20071 ). Such 



models simulate how planetary embryos accrete water-rich km-sized planetesimals coming 
past the snow line. 

Although it is currently not possible to est imate the occurrence rate of hot desert 



worlds, we infer from N-body simulations (e.g., iRaymond et al 



20041 ) that metal poor 



and/or slowly rotating stars could preferentially host such planets. The most important 
factor influencing the water content of small exoplanets is the surface density of solids 
around the snow line. If the surface density is low, dry Earth-sized planets form. Two main 
factors influence the surface density of solids: the total disk mass, and the metallicity. The 
disk mass is influenced by the angular momentum of the protostellar core. If the core rotates 



slowl y, most mass falls directly onto the protostar, and t 



[e.g., 



Yorke et al 



1993 : 



Matsumoto et al.iil997 : 



le disk mass remains moderate 



ZhuetaL 



2010l ). As angular momentum is 



conserved, a slowly rotating core produces a slowly rotating star. Thus the stellar rotation 
period might be indicative of the exoplanetary water content. Metallicity also influences 
the surface density of solids past the snow line. Thus hot desert worlds could preferentially 
form around low metallicity stars. The occurrence rate of low-r nass planets does not seem 



Buchhave et al 



20121 ) ■ but 



to have a strong metallicity dependence (unlike for hot Jupiters 

a statistical sample of characterized exoplanet atmospheres could in the future show us that 

the exoplanetary water or volatile content correlates with the metallicity. 
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4.2. Observables 

Do desert worlds have a hallmark signature in their spectra? Can we identify a 
habitable desert world? The surface conditions (such as surface pressure, temperature, 
and the mixing ratio of atmospheric constituents) must be constrained by observations to 



unambiguously identify a habitable world. We use transmi ssion spectroscopy to i 



lustrate 



the difficulty in this process. We show the effective height (iKaltenegger fc Traubl 120091 ) of 
three types of atmospheres in Fig. [121 our nominal hot desert world with 10 m/s^ surface 
gravity (Tabled]); the effective height of a CO2 dominated Venus-like atmosphere with a 
water mixing ratio of 3 x 10~^; and the effective height of an Earth-like atmosphere with 
N2/CO2/H2O (other gases are removed for easy comparison). The spectra are computed at 
high resolution {X/6X = 10^) and binned down to a resolution of 1000, whic h is comparable 



to th e medium resolution of NIRSPEC on the James Webb Space Telescope (JBagnasco et al. 



20071 ). The illustrated transmission signals are noise- free, but a typical signal -to-noise ratio 



is exp ected to be ~10 for water features at 2 micron for systems out to 12 pc ( JDeming et al. 



20091). 



It might not be possible to identify a habitable exoplanet with certainty. The 
greatest challenge is to constrain the surface pressure and the mixing ratio of water in 
the surface regions. Water vapor is abundant in the atmosphere close to the surface, if a 
surface reservoir is present. However, transmission spectroscopy typically constrains the 
stratospheric mixing ratios only, because the surface region is expected to be optically 
thick at all wavelengths, and thus the transmission signal is, for the most part, unaffected 
by the troposphere. This is a problem, because the surface properties such as the surface 
pressure, and the water reservoir of the planet remain unconstrained by observations. This 
point is illustrated in Fig. [13] because the water features of a hot desert world and an 
Earth-like planet are similar, although the surface relative humidity differs by almost two 
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orders of magnitude. Furthermore, the effective height is above 5 km for both cases, thus 
the surface regions are unconstrained. The largest difference in the transmission signal is at 
wavelengths shorter than 1 micron, because this part of the spectrum is influenced by the 
pressure scale height of the atmosphere. The pressure scale height is 

kT 

H, = -r^, (12) 

where k is the Boltzmann constant, T is the temperature of the atmosphere, and M is the 
mean molecular mass. The mean molecular mass is similar in the two cases, the surface 
gravities are identical. However, the temperatures on a hot desert world are larger than 
on Earth. The average surface temperature of the hot desert world is 370 K in contrast to 
the 290 K surface temperature of Earth. The amplitude of the molecular features differs at 
all wavelengths. This difference is also explained by the large pressure scale height on hot 
desert worlds. 

Although we might not be able to tell with certainty whether a planet is habitable, 
it will be possible to rule out inhabitable planets with confidence. If the stratospheric 
mixing ratios of CO2 and water are retrievable from the transmission spectrum, it is 
possible to distinguish a C02-dominated atmosphere from an N2-dominated atmosphere 



(JBenneke fc Seager 



2012). Variations in the effective height are modest in a CO2 dominated 
atmosphere for two reasons. First, the CO2 absorption lines are wide and the atmosphere 
absorbs significantly even in between band center regions. Second, the scale height of a 
CO2 dominated atmosphere is small because the molecular weight of CO2 is larger than 
the molecular weight of N2. However, if the CO2 and water mixing ratios are small, the 
variations in effective height are large (see Fig. [T3|) . Measuring the relative effective heights 
of band centers and window regions is a way to distinguish C02-dominated atmospheres 
from Earth-like and hot desert world atmospheres. 

We emphasize the previously known result that the transmission signal of low surface 
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gravity planets is easier to detect because the atmosphere has a large scale height (Eq. [T2|) 
and thus the variations in the transmission signal are large. We illustrate the influence of 
surface gravity on the effective height for hot desert worlds in Fig. [HJ 
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Fig. 13. — The effective height of a hot desert world with 10 m/s^ surface gravity (other 
parameters in Table [1]), Venus (with the km height moved to 1 bar), and an Earth- 
like planet (with only N2, water and CO2 in the atmosphere). Medium resolution spectra 
obtained with the JWST will enable us to distinguish a hot desert world from a Venus-like 
atmosphere. The effective height of an Earth-like planet and a hot desert worlds look similar. 
The difference at short wavelength is due to the different scale heights of the atmospheres. 
However, the surface regions remain optically thick at all wavelengths. 

4.3. Model uncertainties 

4.3.1. Clouds 



The effects of clouds on the inner edge distance is presumably small, because we do 
not expect tenuous water clouds covering a large fraction of the surface on dry worlds. 
However, there are some cases where clouds could influence the inner edge distance. Most 
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Fig. 14. — The effective height of hot desert worlds with various surface gravities (for other 
parameters, see Table [1]). The surface gravity influences the pressure scale height of the 
atmosphere. The molecular features are easier to detect, if the exoplanet's surface gravity is 
small, thus the pressure scale height is large. 

notable is clouds forming in the convective equatorial regions of the Hadley cell. Cloud 
formation is initiated by condensation, thus large surface-tropopause temperature difference 
is necessary to form convective clouds. The temperature difference is large, if the CO2 
mixing ratio and/or the surface pressure are large. Even if clouds do form, it is expected 



that clouds form high in the troposphere ( IZsom et al 



20121 ). therefore the clouds might be 



thinner on hot desert worlds than on Earth. Furthermore, it is difficult to a priori estimate 
whether such clouds would reduce or increase the inner edge distance, because the albedo 
effect and the warming efj Fect of high clouds is typically on the same order of magnitude 



(IGoldblatt fc Zahnle 



20111 ). 



Clouds forming on the night side of slow rotators could also influence the inner edge 
distance, if such clouds extend to the day side. If zonal winds shepherd clouds to the day 
side, and the clouds survive to the vicinity of the sub-stellar point, the albedo effect of 
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clouds could push the inner edge closer to the host star. Clouds form predominantly at 
the poles on fast rotators assuming that their rotation axis is close to perpendicular to the 
orbital plane, thus their effect on the inner edge distance is limited due to the large zenith 
angle of incoming stellar radiation. 



4.3.2. ID vs. 3D 

The main assumption in our approach is that the mixing ratio of atmospheric 
constituents and the temperature- pressure profile is described by ID vertical globally 
averaged quantities. If the mixing ratio of a chemical is influenced by processes operating 
on short timescales relative to atmospheric mixing, spatial variations in mixing ratio 
throughout the atmosphere are expected to be large. CO2 i s assumed to be we ll-mixed 



because the carbon-silicate cycle operates on a long timescale ( iWalker et al.iil981i) . However, 
the mixing ratio of H2O could show large horizontal variations because the surface waters 
are expected to be confined to a small portion of the surface on hot desert worlds, and 
precipitation is not anticipated on the day side. Such mixing ratio variations could have 
important implications to planetary climate and the inner edge distance because regions 
with below average relative humidities act as global coolants. Two examples are the dry 
downwelling part of the Hadley cell or the Foehn wind. The optical depth in these columns 
is reduced and more infrared radiation escapes to space thus cooling the planet. Areas 
with above average relative humidities are expected where precipitation occurs (night side 
or poles). The enhanced greenhouse effect in these areas increases the cooling timescale of 
Eq. [7] thus reduce the overall horizontal temperature gradients. Multidimensional coupled 
radiative and circulation models are necessary to study these effects on the climate. 
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4-3.3. Relative humidity profile 

Other crucial assumption is the constant tropospheric relative humidity profile. The 
relative humidity profile determines the column density of water vapor in the troposphere 
and thus the radiative properties of the atmosphere. The relative humidity profile on 
Earth decre ases with altitude being 70% on the surface and a couple of percentage at the 



tropopause ( jManabe fc Wetherald 



19671 ). The relative humidity profile is determined by the 
water cycle: surface evaporation, atmospheric circulation, condensation, and precipitation. 
It is beyond the scope of this paper to study such effects, but we plan to investigate this 
problem in the future. 



5. Summary and Conclusions 

We studied atmospheric configurations that provide habitable surface conditions while 
minimizing the semi-major axis of the exoplanet. We recalled that liquid precipitation 
is a dominant removal process for atmospheric CO2. Rain is also necessary to keep the 
relative humidity low because it acts as a sink for both water vapor and CO2. Therefore, 
we provided a lower limit on relative humidity required for liquid precipitation. We studied 
the atmospheric circulation requirements to initiate condensation and precipitation either 
on the night side or at the poles of exoplanets. The water loss timescale was also estimated 
to assess the life-time of the surface water reservoir. Note that we focused exclusively on 
the inner edge of the Habitable Zone, because the o u ter edge extends well b eyond the 



N2-O2 Habitable Zone limits of 



dominated atmospheres (jPierrehumbert &: Gaidos 



Kasting et al. 



(Il998h: 



Kopparapu et al. 



torn for H2-He 



20 111 ) , or it can extend to infinity, if the 



exoplanet has a sufficient internal heat source (Rogers et al, in prep.). 



Our main results are summarized here: 
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• We estimate that the fundamental inner edge of the Habitable Zone is at 0.5 - 0.6 AU 
around a solar-like star, well within the orbit of Venus for a variety of atmospheric 
scenarios. 

• Planets with dry atmospheres are habitable close to the host star. We estimate that a 
relative humidity of 1% could be sufficient to initiate hquid water precipitation. If the 
relative humidity is smaller than 1%, precipitation might occur predominantly in the 
form of snow that is unable to remove CO2 from the atmosphere if the temperature is 
above the freezing point of CO2. 

• The water loss timescale is longer than 10 billion years, if the atmospheric CO2 level 
is above 10~^ and the surface pressure is more than 1 bar. As we argued in Sect. 12. 3[ 
if plate tectonics does not recycle carbon, the atmosphere and the mantle could be 
depleted in CO2 on a billion year timescale, and carbon remains locked in the form of 
carbonates in shallow water regions. Therefore plate tectonics or other processes are 
necessary to either keep the CO2 level balanced, or to reduce the stratospheric water 
mixing ratio by other means. 

• If the surface pressure on slow rotators and tidally locked planets is on the order of 
1-10 bars, it is plausible that zonal winds enable a cold night side (dew point) and 
precipitation. For large surface pressures (e.g., 100 bars), the heat capacity of the 
atmosphere is large and it prevents strong cooling and temperature differences. 

• Liquid water precipitation is plausible at the poles of fast rotators, because the 
required baroclinic diffusion coefficient to initiate precipitation is large for all 
atmospheric scenarios considered here. 



• The number of potentially habita b 
as compared to 



Kopparapu et al 



e con firmed exoplanets increases by a factor of 2-3 



(I2OI3I ). The potentially habitable exoplanets are 
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Gliese 581 d, e, f, Gliese 667Cc, Gliese 163c, HD 20794d, HD 40307g, tau Ceti d, e, f, 
and Kepler 22b. 

• The surface regions are not constrained by transmission spectroscopy even though the 
water mixing ratio on hot desert worlds is smaller than on Earth. 

The first confirmed habitable exoplanet will likely be a hot desert world - as introduced 
in this study - because it is easier to characterize their atmospheres. Such close-in exoplanets 
orbiting small stars have a large transit probability, and a large transit over rotation period 
ratio that increase the observational window for a given mission. Characterizing low mass 
desert worlds around nearby M dwarfs is especially likely because the atmospheric scale 
height increases for low-mass planets while the planet-to-star area ratio remains favorable. 
Therefore future modeling studies will aim to better understand the atmospheric circulation, 
the water cycle, and geophysical processes of such exoplanets. Observational campaigns 
conducted to characterize Earth-like exoplanet atmospheres should consider all super-Earth 
or Earth-sized exoplanets residing further out their host star than the fundamental inner 
edge limit. 

A. Zsom is grateful for discussions with Nikole Lewis, Renyu Hu, Vlada Stamenkovic, 
Ramses Ramirez, Daniel James Cziczo, Ray Pierrehumbert, James Kasting, and Felipe 
Gerhard. A. Zsom was supported by the German Science Foundation (DFG) under grant ZS 
107/2-1. J.d.W. acknowledges support from the B.A.E.F. (Belgian American Educational 
Foundation) in the form of a fellowship. 
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A. Model validation 



Thermal emission We compare our radiative transfer code against the radiative 
transfer (RT) module of t he Community C 



for Atmospheric Research (JKiehl et al 



imate Model (ccm3) of the National Center 
19981 ) ■ The ccrnS model was developed for climate 
change studies, thus its radiative transfer module is more accurate than our code in the 
temperature, pressure, water and CO2 mixing ratio ranges typically encountered in the 
atmosphere of Earth. Several effects are included in the ccmS RT model that we ignored 
to keep our code simple. Some of these effects are: properly treating overlapping lines of 
different absorbers (we use the assumption of random overlapping), calculating the IR flux 
in subdivided atmospheric layers to resolve the temperature gradient through the layer (we 
assume that the layers are isothermal) . We show by the end of this section that ignoring 
these effects results in a relative error on the order of a few percent in the thermal IR fluxes. 
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Fig. 15. — The outgoing thermal flux of a dry N2-O2 atmosphere with 1 bar surface pres- 
sure, 273 K surface temperature, and various CO2 mixing ratios. The solid line represents 
calculations performed with the Community Climate Model (ccm3) of the National Center 
for Atmospheric Research. Thermal flux values calculated with our model are shown with 
'+' signs. The relative error between our results and the ccmS model is below 2%. 
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Fig. 16. — Left: The outgoing thermal flux of a moist adiabat atmosphere for four levels 
of relative humidity and various surface temperatures without CO2. Right: The outgoing 
thermal flux of moist adiabat atmosphere profiles with 10^ ppmv CO2 added. The solid lines 
correspond to ccm3 calculations, the '+' symbols show our results. The relative error is less 
than 7% in all cases. 

We calculate the outgoing thermal emission {OTE) of all-troposphere atmospheres 
with various temperatures, relative humidities and CO2 mixing ratios using our model 
and ccm3 RT. The first comparison is performed using a dry N2-O2 atmosphere with 273 
K surface temperature and various CO2 mixing ratios. The results are shown in Fig. [T5l 
The solid line represents the ccmS RT results, the '+' signs correspond to our calculations. 
Although our fluxes are typically larger than the ccm3 fluxes, the relative error between our 
results and ccm3 is less than 2%, thus our model agrees well with ccmS. 

Next, the OTE for an N2-O2 all-troposphere atmosphere is calculated with four 
different values of relative humidity: $ = 0, 10, 50, and 100 %, but without CO2 (see Fig. 
[T6k). If the relative humidity is zero, there is no greenhouse gas in the atmosphere, thus the 
corresponding fluxes are generated by the black body emission of the surface. At surface 
temperature below 240 K, the water content of the atmosphere is so low that only a small 
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amount of the IR flux is held back. Fig. [T6k shows how increasing the relative humidity 
lowers the outgoing thermal emission at the top of the atmosphere. The maximum relative 
error observed between our and the ccmS fluxes is 7%, which we find acceptable. 

The next level of complexity is to include both water vapor and CO2 in the atmosphere. 
Fig. [TBb shows the OTE of the atmosphere profiles used in the previous paragraph with 
1% CO2 added. The outgoing thermal flux is significantly reduced by CO2. As before, the 
maximum relative error is an acceptable 7%. 
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